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Abstract 

Regular star formation is thought to be inhibited close to the massive black hole (MBH) in the 
Galactic center. Nevertheless, tens of young main sequence B stars have been observed in an isotropic 
distribution close to it. These stars are observed to have an apparently continuous distribution from 
very close to the MBH < 0.01 pc and up to at least ~ 0.5 pc, suggesting a common origin. Various 
models have been suggested for the formation of the B-stars closest to the MBH (< 0.05 pc; the S- 
stars) . typically involving the migration of these stars from their original birthplace to their currently 
observed position. Here we explore the orbital phase space distribution of the B-stars throughout the 
central pc expected from the various suggested models for the origin of the B-stars. We find that 
most of these models have difficulties in explaining, by themselves, both the population of the S-stars 
(< 0.05 pc), and the population of the young B-stars further away (up to 0.5 pc). Most models grossly 
over-predict the number of B-stars up to 0.5 pc, given the observed number of S-stars. Such models 
include the intermediate-mass black hole assisted cluster inspiral scenario, Kozai-like perturbations by 
two disks, spiral density waves migration in a gaseous disk, and some of the eccentric disk instability 
models. We focus on one of the other models, the massive perturber induced binary disruption, which is 
consistent with both the S-stars and the extended population of B-stars further away. For this model 
we use analytical arguments and N-body simulations to provide further observational predictions. 
These could be compared with future observations to further support this model, constrain it or 
refute it. These predictions include the radial distribution of the young B-stars, their eccentricity 
distribution and its dependence on distance from the MBH (higher eccentricities at larger distances 
from the MBH) , as well as less specific expectations regarding their mass function. 

Subject headings: 
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1. INTRODUCTION 

In recent years high resolution observations have re- 
vealed the existence of many young O and B stars 
in the galactic center (GC). Accurate measurements of 
the orbital parameters of these stars give strong evi- 
dence for the existence of a massive b lack hole (MBH) 
which governs the dyna i iiics in the GC (iGhez et al]ll998l : 
lEisenhauer et all l2005t iGillessen et al.ll2009D . Most of 
the young stars are observed in the central 0.5 pc around 
the MBH. The young stars population in the central 
~arcsecond (0.05 pc), the so called 'S-stars', contain only 
young B-stars (likely masses of 7 — 15 M©), in apparently 
isotropic distribution around the MBH, with relatively 
high eccentricities (consiste nt with a thermal di stribu- 
tion, with 0.3 < e < 0.95; IGillessen et"an 120091) . The 
young stars outside this region contain many 0-stars 
and present markedly different orbital properties, with 
a top heavy mass function and a disk (or disks) like con- 
figuration(s) jLevin fc Beloborodo-(l l2003l;lLu etal]|2009l; 
iBartko et all 120091 l20ldi r and much more circular or- 



bits. The sharp difference and discontinuity between 
the 0-stars population outside the central 0.05 pc, and 
the S-stars population inside this region strongly sug- 
gest a different ori gin for these two populations (see also 
IBartko et aLll2010[ ). 

Recently, it was found that the isotropically distributed 
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B-stars population extends (apparently) continuously be- 
yond the c entral arcsecond and up to at least 0.5 pc from 
the MBH (jBartko et al.ll2010t ). Although many studies 
have explored the origin of the 'S-stars' (see Perets et al] 
120091) ■ the observations of the extended population may 
shed new light on the origin of these stars. Given the 
seemingly continuous nature of the B-star population, 
it is likely that this population of both the S-stars and 
the extended population at larger distances (or at least 
most of it), has a single origin. In the following we make 
this assumption and study possible scenarios for a sin- 
gle common origin of all the GC B-stars and their later 
evolution. We discuss and suggest specific observational 
constraints on the various models. We then mainly focus 
on one of these m odels, the massive perturbers induced 
binary disruption ([Perets et al.ll2007[ ). and use analytical 
arguments and N-body simulations to provide specific 
predictions and observational signatures expected from 
this model. 

2. ORIGINS AND EVOLUTION 

Regular star formation in the close region near the 
MBH is inhibited by the tidal forces which destroy any 
potential giant molecular cloud progenitor. For this rea- 
son, young stars in this region are thought to either form 
elsewhere, further away from the MBH, and then mi- 
grate rapidly closer to the MBH; or form in-situ near 
the MBH, in an irregular form of star formation, most 
likely through instabilities in a gaseous disk (see below). 
After their formation/capture, stars in the GC environ- 
ment can dynamically evolve. Each of the scenarios for 
the B-stars origin produces a different and unique distri- 
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bution of their initial orbital configurations. These con- 
figurations, however, can later dynamically evolve due 
to relaxation processes. The currently observed prop- 
erties of these stars therefore reflect the combination of 
the initial configurations and their later evolution. In the 
following we discuss the various models for the origin of 
the B-stars. We briefly review each of them, discuss their 
possible caveats, and detail dynamical constraints which 
could be used to exclude or support these models. As we 
shall see, simple efficiency arguments show that most of 
the suggested models for the B-stars origin can already 
be strongly constrained, if not excluded, by current ob- 
servations. 

Note that although we focus on the origin of the B-stars 
in the GC, in the following we also mention scenarios for 
the disk O and WR-stars origin, since their population 
is closely related to that of the B-stars in some of the 
models we discuss. 

2.1. In-situ star- formation from inf ailing gaseous 
clumps 

It was suggested that the young OB stars in the GC 
were formed in-situ in the central pc a few million years 
ago in gaseous disks and/or strea ms formed fro m infalling 
and/or colliding gaseous clumps (Morris 1993). Analytic 
calculations and simulations ( Navaks hin fc Cuadra jOOS,; 
Levin 2007; BonneU & Rice 2008: IWardle fc Yusef-Zadehl 
2008i : iHobbs fc Navakshin ,2009i) have shown that stars 
could form in such fragmenting clumps to produce stel- 
lar disks and/or other coherent structures, in the region 
of a few 0.01 pes up to a few 0.1 pes from the MBH. These 
could, in principle, be consistent with the observed young 
stellar structures (clockwise disk and a secondary in- 
clined structure) . This scenario could also explain the ra- 
dial distribution of stars in the observed stellar disk and 
possibly their eccentricities. Moreover, some of the pro- 
posed models produce a top heavy mass function (MF) 
for the ne wly formed stars, po ssibly consistent with ob- 
servations ijBartko et al.l[2010l ). None of these, however, 
produces the population of B-stars within 0.05 pc of the 
MBH. 

One should note that our poor understanding of the 
initial conditions and the star formation processes in the 
GC, allows for a wide parameter space to choose from, 
which, naturally, raises difficulties in constraining (or 
strictly falsifying) such models. Nevertheless, the robust- 
ness of producing some sort of star formation in the GC 
region under a variety of conditions explored in the liter- 
ature suggest these models as the currently most promis- 
ing scenarios for the origin of the young stellar disk and 
stellar structures (although less likely the origin of the 
isotropic B-stars population, including the S-stars, as dis- 
cussed below). 

2.2. In-situ star-formation followed by rapid migration 

None of the models for in-situ star formation in the 
central pc suggests the formation of stars as close to the 
MBH as the S-stars, or the apparent existence of two dis- 
tinct young O and B stellar populations. Producing both 
the disk and isotropic B-stars population in the same sce- 
nario (and in particular the inner S-stars) requires some 
fine-tuning. Specifically, comparing the top heavy mass 
function of the disk stars, as suggested by current ob- 
servations, with the more regular mass function of the 



isotropic B-stars disfavors a common origin. One there- 
fore requires a selective process which works differentially 
on stars of different masses. In any case, an additional 
process would be required for the migration of stars from 
the outer to inner regions in the central pc in order to 
produce the S-stars close to the MBH. 

Two-body relaxation processes work, in principle, dif- 
ferentially on stars of different masses, through mass seg- 
regation (energy equipartition) . For example, relaxation 
in an isolated stellar disk would somew hat segregate the 
more massive stars into a thinner disk ([Alexander et al.l 
I2007t iPerets et al] l2008bt iLockmann et al l 120091) . pro- 
ducing mass stratified populations. However, this effect 
is relatively small. Moreover, the most important com- 
ponent for relaxation in the GC is likely to be the stellar 
black holes population in the stellar cusp dPerets et al.l 
I2008bt ILockmann et al.ll2009t iPerets et all I2009D. which 
is dominated by resonant relaxation ( Perets et al.ll20091) . 
Segregation into two different stellar populations of dif- 
ferent mass regimes is not likely to ensue in this case. 
Moreover, such scattering of stars could not produce the 
popula tion of S-stars closer to the MBH (jPerets et al.l 
l2008bt also Perets et al., in prep . ). Encounters of bin ary- 
single stars (|Cuadra et al.ll2008t IPerets et alll2008aD can 
also only have a small effect in producing the isotropic B- 
stars population from a thin stellar disk. P erturbations 
by massive perturbers (jPerets et al.l l2007f ) su ch as in- 
falling intermediate - mass bl ack ho les (IMBHs) (|Yu et al.l 
|2007; Gualandris fc Merrittll2009l ) or other stellar disks 
((Lockmann et al.ii2009l ; also Gualandris et al., in prep.) 
do not differentiate between stars of different masses ei- 
ther. 

2.2.1. Efficiency constraints on in-situ formation of the 
B-stars in the central pc 

Given the lacking suggestion for a mass differentiat- 
ing process, a different route can be taken to explain the 
different stellar populations. We can suggest that two 
distinct epochs of in-situ star-formation occurred. Two 
such epochs would still require migration of some of the 
stars from the disk to the inner region, where stars do 
not form in-situ. One should mention, in this context, 
that in some cases a different mass function was found for 
stars formed in different struct ures in the same simula- 
tion C Hobbs fc Navakshinll200a . We can therefore either 
suggest two epochs of star formation happening at differ- 
ent times, or a single epoch producing two distinct popu- 
lations. In either case, the rapid migration producing the 
S-stars should affect only one of the stellar populations 
formed. 

Two distinct star formation epochs could naturally pro- 
duce more massive and less massive stellar populations. 
Even if both epochs produced stellar populations with 
the same initial mass function, the most massive stars 
from the first epoch might already have ended their life, 
leaving behind a stellar population of less massive stars. 
This scenario could be constrained by observations. We 
may suggest an interesting and very general requirement 
applicable to any scenario in which the S-stars formed 
far from their current positions. In any such scenario 
only a fraction fmig of the stars formed in some external 
region (e.g. the stellar disk, the central pc, or regions 
outside the central pc) finally migrate to become S-stars. 
For a given number of observed S-stars, Ns, the parent 
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external population should be 1/ fmig larger to have 

-^par s I fmig (1) 

stars (compa re with similar constraints derived by 
iPeretj l2009a|) . Cases in which f^ig is small could be 
strongly constrained by such requirement. For exam- 
le, the 'billiard' model for the origin of the S-stars 
Alexander fc Livioll2004D in which stars from the cen- 
tral pc are captured close to the MBH through ex- 
changes with SBHs close to the MBH, is disfavored 
since the population of similar B-stars in the central 
pc is too s mall to accommodate the required parent 
population ([Paumard eTall 12001 . Similarly, one can 
turn to models of a disk origin for the S-stars, such 
as the Kozai-like perturbatfo ns in the two disks sce- 
nario (jLock mann et al.ll2009l). the eccentric disk insta- 
bility model (Madigan et al. 200^, or the spiral density 
wave model (Griv 2010). All of these models, irrespec- 
tive of details, suggest that the S-stars formed as part of 
a stellar disk, up to a few 0.1 pc away from the MBH. 
A small fraction of the stars formed in the disk migrated 
through some process to later on become the currently 
observed S-stars. The current number of B-stars inferred 
in the central 0.5 pc, but outside the central 0.05 pc 
where the S-stars reside, is comparable, and likely some- 
what smaller than the number of S-stars. These models 
are therefore required to have a very high migration ef- 
ficiency in order to explain the origin of the S-stars. In 
table [1] we list the expected number of B-stars in the cen- 
tral 0.05 pc and those up to ^ 0.5 pc, according to the 
various models. We compare these numbers with the ob- 
servationally inferred numbers. We also list some of the 
non in-situ formation models, which are discussed later 
on. Most of the current formulations of the models are 
likely excluded, as can be seen from the table. Whether 
similar models could still be consistent with observations 
under somewhat different conditions than currently stud- 
ied, would require additional investigation. 

2.3. External star-formation followed by rapid migration 

Another set of scenarios suggest that the GC young 
stars formed outside the central pc, where conditions are 
less hostile to regular star formation. These scenarios 
could be used to explain either the origin of the stellar 
disk, the origin of the isotropic B-stars, in particular the 
S-stars, or both. Note that the extension of S-stars dis- 
tribution beyond the central 0.05 pc is only a recent ob- 
servational development, nevertheless, some of the sug- 
gested models discuss the possible existence of such an 
extended distribution. 

2.3.1. Cluster and IMBH infall 

An infall of a young stellar cluster (with or with- 
out an IMBH) into the central pc was suggested 
as an alternative scena rio fo r the origin of the GC 
youn g stars ( 'Gerhard 200lt_lH ansen k, Milosavlievid 
20031: iK im fc Morris 2003; Porteeies Zwart et al.l l2003l; 
Kim et al. 2004: Levin et al. 2005; Giirkan & Rasio 2005; 
Berukoff fc Hansen. ,2006: ,Fuiii et al.,,2008, . ,2009il . Such 
dissolving cluster is likely to form a stellar disk-like 
structure, possibly with additional outlying structures 
and/or isolated stars outside the main disk, as ob- 
served in the GC. It would also produce a bias towards 



more massive stars being concentrated in the central re- 
gion of the GC. These massive stars, which were orig- 
inally segregated in the innermost regions of the clus- 
ter, would be the last to evaporate from the cluster, 
i.e. in the central most regions closer to the MBH. In 
addition, an IMBH infall has the potential to deplete 
the i nner regions of the GC, and form a core-hke st ruc- 
ture (jLockmann fc Baumgar"d^ '20081: jMerrittI I2009D . as 
possibly observed (iBuchholz et al. 20"09|: iDo et al.ll200a 
iBartko et a"ni2010[ ). lYu et al.l (|2007l) also invoke the ex- 
istence of an IMBH with mass > lO'' Mq to help pro- 
duce higher eccentricities and inclinations for the disk 
stars. An infall of an IMBH could also rapidly isotropize 
the distribution of stars closest to the MBH, possibly 
helping explain the currently obse rved distribut i on of 
the S-stars in the central 0.05 pc (jMerritt et all 120091 : 
iGualandris fc Merrittl[2009l ). These properties make this 
scenario quite attractive. Nevertheless, the inspiral sce- 
nario also has many difficulties, especially in regard to 
the origin of the B-stars on which we focus here. We 
shortly discuss these in the following. 

Infalling clusters may not be able to inspiral in the 
appropriate time win dow for producing the stellar disk. 
iKim fc Morris! (|2003D show that a cluster (without an 
IMBH) would need to be extremely dense or b orn within 
the ce ntral 5 pc to inspiral fast enough. iKim et aP 
(|2004[ ) find that the presence of an IMBH lessens the 
requirements on the cluster density , if it is as massive as 
10% of the cluster mass, however iGtirkan et al.l ()2004l ) 
find that typical masses for IMBHs formed in simula- 
tions of dense clusters (without stellar evolution) are 
of the order of 10"'^ of the cluster mass. The forma- 
tion of appropriate IMBHs in star clusters may there- 
fore be very difficult. Moreover, even the formation 
of lower mass IMBHs through collisions in dense clus- 
ters is found to be difficult when stell ar e voluti on and 
wind mass loss are taken into account dY ungelson et al] 
l2008l:lGlebbeek et al.ll2009l : IVanbeveren et al...200a ). The 
IMBH cluster infall model is therefore unfavorable given 
our current understanding. The mass and orbital sepa- 
ration of a possible IMBH in the GC are constrained by 
radio observations of SgrA* (see IGualandris fc Merrittj 
l2009l for a recent overview of these constraints). 

Another potential issue with the infall model is that 
an inspiraling cluster is likely to leave most of its stars 
behind as it dissolves, whereas very few young stars 
are observed outside the central 0.5 pc. Simulations 
of an infalling cluster hosting an IMBH ( Levin et al.l 
l2005t iBerukoff fc HansenI l200a iFuiii et all [2010D show 
that the young stars are stripped from the cluster be- 
fore the IMBH reaches the central 0.05 pc. The young 
stars closest to the MBH (the S-stars) are therefore not 
likely to directly originate from such a scenario. Starting 
from a cluster of 50 stars strongly bo und to an inspiraling 
IMBH (separation of up to 0.06 pc) . IBerukoff fc HansenI 
([2006) showed that at most 1-2 stars could reach dis- 
tances of 0.04 pc from the MBH. In other words, 96 — 98 
percent of the stars were stripped off when the IMBH 
was at distances of a few 0.1 pc from the central MBH. 
In terms of the "migration efficiency" defined above, we 
get fmig < 0.04, where we have replaced the parent pop- 
ulation of the newly formed stars in the disk with the 
parent population of the stars originally clo sely bound 
to the IM13H. Similar results are obtained bv lFujii et al.l 
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Model 


fmig 


Observations+ 


Expected * 


Observations+ 


Ref. 






0.05 pc 


< 0.5 pc 


< 0.5 pc 




Massive perturbers (GMCl) 


~ 0.25 


66 ± 15 


~ 264 


212 ±50 


This paper 


External binaries (Empty loss cone) 


~ 0.76 


66 ± 15 


~ 87 


212 ±50 


This paper 


External binaries (Full loss cone) 


~ 0.12 


66 ± 15 


~ 539 


212 ±50 


This paper 


Eccentric disk instability {e-mit > 0.7) 


> 0.16 - 0.321 


66 ± 15 


< 206 - 412 


212 ±50 


Madigan et al. 2009 


Eccentric disk instability {cinit =0.6) 


~O.O7-O.i4i 


66 ± 15 


~ 472 - 943 


212 ±50 


Madigan ct al. 2009 


Eccentric disk instability'^ {einit < 0.5) 


< 0.015 - O.O3I 


66 ± 15 


> 2200 - 4400 


212 ±50 


Madigan ct al. 2009 


Two disks Kozai"^ 


< 0.045 


66 ± 15 


> 1467 


212 ±50 


Lockmann ct al. 2009 


Cluster + IMBH inspiral 


< 0.04 


66 ± 15 


> 1650 


212 ±50 


IBerukoff & Hansen 2006: Fuiii et al. 2010| 


Spiral density wave* 


~ 0.01 


66 ± 15 


~ 6600 


212 ± 50 


lGrivil20lOI 



^Inferred from observations (stars with K magnitude 13 < 14 < 16.5). [Sartko et al.l i2010D also H. Bartko, private communication, 2010 
*The expected number normalized to get ~ 53 stars at < 0.05 pc 

^The range corresponds to a range of total binary fraction in the disk, taken here to be 0.5 — 1. 
^Several models studied, most efficient studied model is shown (in this case a disk with einu — 0.5) 



"^Several models studied, most efficient studied model is shown (4c in lLockmann et al.ll2009l) 
''The total number of S-stars goes like (Aa/ Act = (0.05/0.5)'^ ~ 0.01, seelGrivll2010l 



(|201Cl( ). but under somewhat different conditions, who 
find one captured S-star (<0.05 pc) for 21 stars stripped 
at < 0.4 pc (M. Fujii, private communication, 2010). For 
the current population of S-stars to originate from such 
a scenario one requires (Eq. [ij Np = Ng/ fmig ^ 1000 
B-stars to be stripped from the cluster and currently ex- 
ist in the central 0.5 pc; this is much larger than the 
currently inferred number. 

It is also difficult to envision a super-massive cluster 
containing > 1000 massive B-stars within 0.06 pc of the 
IMBH for this scenario to have the necessary apriori con- 
ditions. 

Even though none of the existing simulations of IMBH 
inspiral have been able to pr oduce the S-clust er, includ- 
ing the most recent ones by iFujii et alJ (|2010l ) in which 
massive stars are carried down to distances no smaller 
than a fe w 0.1 pc f r om th e MB H, we do note tha t only 
[Berukoff fc HansenI (|2006[ ) and iFuiii et al.l ()2010[ ) have 
studied in depth the IMBH infall scenario in the con- 
text of the S-stars. Future studies may explore different 
conditions than those used there and may bring new in- 
sights on this subject. Similar efficiency arguments as 
discussed here could then be applied to test the consis- 
tency of such simulations with current observations. 

2.3.2. Binary disruption 

Another scenario where young stars are formed far 
from the MBH and then migrate close to it is the bi- 
nary disruption scenario. A close pass of a binary star 
near a MBH results in an exchange interaction, in which 
one star is ejected at high velocity, while its co mpan- 
ion is captured by the MBH and is left bound to it (jHillsl 
Il988f l. Such interaction occurs because of the tidal forces 
exerted by the MBH on the binary components. A young 
binary star could therefore be formed outside the central 
region and later be scattered onto the MBH on a highly 
radial orbit leadi ng to its disrupt i on. S uch a scenario 
was suggested byl Gould fc QuillenI (|2003[ ) to explain the 
origin of the star S2. In order for the capture rate of such 
stars to explain the current observation of all the GC B- 
stars, rapid relaxation processes are required for the bi- 
naries to be scattered onto the MBH. Such a model, sug- 



gested bv iPerets et al.l ()2007() . which takes into account 
scattering by massive perturbers outside the central 1.5 
pc (such as giant molecular clouds and clu mps observed 
in th e GC region and other gala ctic nuclei; iPerets et alJ 
l2007t IPerets fc Alexanded 120081 ) could possibly account 
for the observed number of B-stars. 

The binary disruption scenario leaves the captured 
stars on highly eccentric orbits (e > 0.95), and further 
dynamical evolution is required in order to explain their 
currently observed more relaxed eccentricity distribution. 
Study of their evolu tion, which is driven by resonant 
relaxation processes (jRauch fc Tr emaind ll996t ). suggest 
that indeed the more relaxed, almost thernial ecc entric- 
ity distribution of the S-stars (jPerets et al.ll200"9l) could 
be consistent with their evolution from a much higher 
initial eccentricity. 

Note that this scenario, like the disk origin models for 
the S-stars, discussed above, can be constrained by ob- 
servations of the parent population from which the S- 
stars or iginate. Current observations do not exclude this 
model (jPerets et al.l 120071 ; lPerets^l2009a^ . Future obser- 
vations searching for young stars in these regions should, 
in principle, give better constraints. Nevertheless, other 
observational signatures may be more easily verified or 
refuted. We discuss these in the following section which 
focuses on this scenario. 

3. EVOLUTION OF CAPTURED STARS FROM BINARY 
PROGENITORS OUTSIDE THE CENTRAL PC 

In the previous section we discussed various models for 
the origin of the B-stars. As discussed above, most of the 
suggested models can be strongly constrained and even 
excluded by simple arguments based on current observa- 
tions. In the following we focus on the binary disruption 
model, not yet excluded by sim ilar arguments. In a pre- 
vious study (jPerets et al.ll2009h we focused on the S-stars 
origin, and used analytic arguments and TV-body simula- 
tions to show that a binary disruption origin followed by 
the (resonant) relaxation of the captured stars is consis- 
tent with their currently observed dynamical properties. 
Here we extend this study to explore the dynamics of the 
isotropic population of B-stars throughout the central 
pc. We first describe in more details the binary disrup- 
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tion scenario. We then discuss the radial distribution of 
the captured B-stars according to the model and how it 
compares with observations. We also shortly discuss the 
mass function of such stars. Finally, we use A'^-body sim- 
ulations and analytic arguments to find the expected ec- 
centricity distribution of the B-stars and its dependence 
on the distance from the MBH. We show that the latter 
distribution is likely to have quite unique properties and 
could therefore serve as a strong constraint on the binary 
disruption model. 

3.1. The binary disruption scenario: details 

Young stars forming far from the MBH could be part 
of binaries. If such a binary is scattered into a highly 
eccentric orbit, it can pass the MBH very closely. A close 
pass of a such binary near the massive black hole may 
result in an exchange interaction, in which one star is 
ejected at high velocity, while its companion is captured 
by the MBH and is left bound to it in a tight orbit. Such 
interaction occurs because of the tidal forces exerted by 
the MBH on the binary components. Typically, a binary 
(with mass, M;,, and semi-major axis, ah), is disrupted 
when it crosses the tidal radius of the MBH (with mass 
M,), given by 



n 



(2) 



The capture probability and the semi-major axis dis- 
tribution of captured stars were estimated by means of 
simulations, which show that most binaries approaching 
the MBH within the tidal radius rtja) (Eg. [H) are ac - 
tually disrupted (jHillj 119911 . 1 19921 ; iBromlev et al.ll2006D . 
The harmonic mean semi-major axis fo r three-body ex- 
changes in equal mass binaries is ()Hill£lll991i) 
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0.56 



M. 



2/3 



0.56 



M. 

Afhin 



1/3 



ru (3) 



where a is the semi-major axis of the infalling binary and 
oi that of the resulting MBH-star "binary". Most values 
of Oi fall within a factor 2 of the mean. This relation 
maps the semi-major axis distribution of the infalling 
binaries to that of the captured stars: the harder (closer) 
the binaries, the more tightly bound the captured stars. 

3.2. Radial distribution of the B-stars 

As discussed above, the initial semi-major axis of the 
captured star around the MBH is linearly related to its 
binary progenitor separation (Hills 1991). The binary 
disruption therefore maps the the distribution of the bi- 
naries separations into the radial distribution of capture 
orbits. However, the binary disruption rates could also 
be dependent on the binary separation, implying that the 
cross section for disruption and hence the size of the loss 
cone increases with separation. Nevertheless, the latter 
is true only for the so-called full loss cone regime (for 
a detailed discussion on disru ption rates and loss cone 
regimes see IPerets et al.]|2007D . 

In the empty loss cone regime, the disruption rate F is 
independent of the size of the loss cone, and only depends 
on the relaxation time scale of the system 

dV dNun{a)/da 

—— oc 

da tr 



(following Eq. 17 in lPerets et al.ll2007l) . where Nbm{a) is 
the number of binaries with separation a and tr is the re- 
laxation time of the system. Since the separation of the 
captured star is linearly related to the progenitor stel- 
lar binary separation, the radial distribution of captured 
stars follows the binary separation distribution. 
In the case of the full loss cone regime, 

dT dNun{a) dNun{a) 

— oc q cx a, 

da da da 

(following Eq. 18 in IPerets et al.l [20071) . where q is the 
the tidal disruption radius of a given binary, which is lin- 
early related to the binary separation (see Eq. |5| above) . 
In this case the radial distribution of the captured stars 
does not directly follow the binary semi-major axis dis- 
tribution. 

We therefore find that the two extreme cases of empty 
and full loss cone regimes produce a different disruption 
rate dependence on the binary separation and therefore 
imply different radial distributions for the captured stars. 

The distribution of semi-major axes of massive binaries 
in the Solar-neighborhood follows a log-constant distri- 
bution. In the empty loss cone regime the disruption 
rate is independent of the binary separations. There- 
fore, the radial distribution of captured stars in the GC 
should directly correspond to the binary separation dis- 
tribution, according to Eq. |3l We therefore expect a log- 
constant, or a cumulative distribution of iV(< r) cx logr 
for the captured B-stars, under the assumption that the 
GC binaries have similar properties to those in the Solar 
neighborhood. In the full loss cone regime, the binary 
disruption rate is linearly dependent on the binary sep- 
aration, and wider binaries have a higher chance of be- 
ing disrupted compared to shorter period binaries. We 
therefore expect a cumulative distribution which goes like 
A^(< r) oc r. 

In fig. [l]we show the cumulative distribution of semi- 
major axes for captured stars up to 0.5 pc, normalized 
to have ^ 40 B-stars within 0.05 pc, fai rly consistent 
with observations (which infer 66 ±15 stars: [Bartko et aTl 
(2010); also H. Bartko, private communication, 2010). 
Also shown in fig. |T]is the expected distribution of cap- 
tured B-stars in the scenari o of binary disrup tion in- 
duced by massive perturbers (jPerets et al.ll2007h . which 
produces ~ 40 B-stars up to 0.05 pc. The distribu- 
tions are obtained usi n g the detailed calculations de- 
scribed in IPerets et al.l ( 20071) . which take into account 
the relevant relaxation processes. The GMCl model de- 
scribed there is assumed for the massive perturber sce- 
nario. A log-constant binary semi-major axis distribu- 
tion is assumed for the progenitor binary population 
(compare with th e log-normal distribution assumed in 
IPerets et al.l |2007|. based on the data available at the 
time; iKobulnickv fc Frveil (|200l . which was later re- 
vised to be a log-constant distribution). This latter sce- 
nario produces a distribution intermediate between the 
extreme empty loss cone regime and the extreme full loss 
cone regime. Interestingly, the captured stars radial dis- 
tribution could therefore give a handle on the type of 
relaxation processes in the GC, given the assumptions 
on the binaries distribution. According to GMCl model, 
the expected number of B-stars in the central 0.5 pc is 
~ 260 stars (when noramlized to 66 at < 0.05 pc). The 
empty and full loss cone regimes, when normalized to 66 
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Fig. 1. — The cumulative distribution of semi-major axes for 
captured stars up to 0.5 pc. Several models are shown, includ- 
ing an empty loss cone, a massive perturbers refilled loss cone 
UPerets et al.ll2007l) . and a full loss cone (see text). The empty 
and full loss cone models are normalized normalized to have ^ 40 
B-stars within 0.05 pc. The massive perturbers model i s not nor- 
malized, and uses the GMCl mode llPerets et ahl [20071 ) assump- 
tions, but with revised binary separtions distribution (see text). 



captured stars at 0.05 pc, predict, respectively, ~ 90 and 
~ 540 stars captured up to 0.5 pc. In order to get a 
simple comparison with other models for the B-stars ori- 
gin, we can define an effective "migration efficiency" for 
this model. In this case we replace the parent population 
with the population of stars captured up to 0.5 from the 
MBH, i.e. strictly speaking, this is not the parent popu- 
lation but the stellar population which one would like to 
probe when comparing the models studied here. We find 



in this case /„ 



0.25 and ^ 0.12 for the empty and 



full loss cone regimes, respectively (see also table [I}. 

3.3. Eccentricity distribution of the B-stars 

The periapse of the captured star in the binary dis- 
ruption sc enario is at r^, and therefore i ts eccentricity is 
very high (|Hillslll991l -. iMiller et al.ll2005h . 

e-1 - rt/ai~ 1 - 1.8{A'hijM,y/^ >0.97 (4) 

for values typical of binaries and the MBH in the GC 
{Mmu = 20 Mq; M. = 4 x 10^ A/©). 

For comparison, the initial eccentricity distribution of 
stars formed in a stellar disk is likely to be biased to much 
more circular orbits. For example, the currently observed 
eccentricity distribution of the young O and WR stars 
observed in the clockwise disk peaks at e ^ 0.35. 

It is interesting to note, however, that the resonance 
relaxation times cales in the GC increase with distance 
from the MBH (jHopman fc Alexandeij [2006af ) . There- 
fore, stars captured/formed further away from the MBH 
are likely to have a less relaxed eccentricity distribution, 
i.e. over time we may expect a correlated eccentricity- 
distance distribution to be produced, with stars more dis- 
tant from the MBH having an eccentricity distribution 
which much more closely resembles the original one. Far 
from the MBH, where the resonant relaxation timescales 
are much larger than the l i fetimes of the B-stars (see 
e.g. iHopman fc Alexandeij l2006af ). stars should keep 
their original eccentricity distribution, i.e. highly eccen- 
tric orbits for captured stars, and likely low eccentricity 



orbits for stars formed in a stellar disk. Stars captured 
a few 0.1 pc away from the MBH, for example, are likely 
to still conserve their initial eccentricities. Closer to the 
MBH captured stars could have a relaxed (thermal ec- 
centricity) distribution even after short times, although 
the captured stars population relaxes to a thermal ec- 
centricity distribution much faster than a st ellar popula- 
tion with initially low typical eccentricity ijPerets et al.l 
MM)). 

3.3.1. Results of N-body simulations 

In order to study the eccentricity distribution of the 
captured stars and its radial dependence more quan- 
titatively we performed detailed iV-body simulations. 
We used the (j> GRAPE direct summation software 
running on gravitySimulator, the 32-node cluster at 
the Rochester Institute of Technology that incorpo- 
rates GRAPE acce lerator boards in each of the nodes 
(iHarfst et al.|[200l . We assume that the main contri- 
bution to the dynamical relaxation of stars in the GC 
comes from stellar black holes (SBHs) . We consider a dis- 
tribution of SBHs following the results of Fokker-Plank 
calculations of t he distribution of 10 Mp SBHs around a 
MBH found by (jHopman fc Alexandeij l2006b( ). We use 
16000 SBHs of 10 M0 with an isotropic radial, r~^, dis- 
tribution between 0.04 pc and 0.8 pc around a MBH of 
3.6 X 10^ Mq (the simulation of stars closer to the MBH 
requires much longer simu lation time; we di scussed the 
evolution of such stars in iPerets et alj |2009() . In order 
not to be affected by boundary effects we focus on stars 

in the range 0.08 0.5 pc from the MBH. In all of our 

simulations we force the MBH to be at rest at the ori- 
gin of the system since the Bro wnian motion is expected 
to be quite low (|Merrittj 120051 ). Note that the fuU cusp 
of low mass stars was not included in our simulations 
(this is beyond our current computatational capabilites). 
Such stars could somewhat accelarate the resonant relax- 
ation processes as more stars contribute to the process. 
However, the larger mass contributing to the cusp poten- 
tial would also produce precession of the stellar orbits, 
thereby producing enhanced quenching of resonant reso- 
nant relaxation. The neglection of the full cusp potential 
therefore suggests that our N-body results likely give an 
upper limit to the relaxation rate due to resonant re- 
laxation. The signature of the initial capture conditions 
might therefore even better preserved than suggested by 
the current simula tions. As in our previous simulations 
(jPerets et al.ll2009| ) general relativisitic precession due to 
the MBH is not included. The latter component, how- 
ever, is not important for stars at the distances from the 
MBH discussed here. 

The "captured stars" followed by our simulations are 
those stars with initial eccentricities of 0.95 < e < 0.99. 

In fig. [2]we show the temporal evolution of the median 
eccentricity of the captured B-stars in the central 0.5 pc, 
at various distances from the MBH. The MS lifetimes of 
7—15 Mq stars are '^10 — 50 Myrs. in the future. 

In figs. [3K-l3t we show the cumulative eccentricity dis- 
tribution in three different distance regimes and at 20, 40 
and 60 Myrs of evolution. Also shown for comparison is a 
thermal eccentricity distribution (/(e) — 2e : cumulative 
c/(e) = e2). 

In fig. 21 we show the median eccentricity of the cap- 
tured B-stars vs. their separation from the MBH at 20, 
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Fig. 2. — The temporal evolution of the median eccentricity 
of the captured B-stars in the central 0.5 pc, at various distances 
from the MBH. 



40 and 60 Myrs of evolution. As can be seen in this 
figure, the binary capture scenario provides a qualita- 
tively unique signature, in which the typical eccentricity 
of the stars rises with distance. The exact relaxation 
rate in the GC is dependent on the number density of 
stars and SBHs in this environment, which is still poorly 
constrained. For our assumptions on the number den- 
sity of SBHs, captured stars a few 0.1 pc away from the 
MBH have very high (> 0.9) eccentricities, whereas stars 
are more relaxed at closer to the MBH, as also found in 
previous stu dies on the capture d population even closer 
to the MBH ( Perets et al.l l200^. Faster relaxation rates 
(e.g produced by a higher density of SBHs and/or the ad- 
ditional contribution to relaxation process by low mass 
stars, not included here) would still produce the quali- 
tative signature of a typical eccentricity rising with dis- 
tance, but may quantitatively show lower eccentricities, 
and a relaxed population up to larger distances from the 
MBH. Conversely, a stellar population with initially low 
(< 0.5) eccentricities, as might be expected from a stel- 
lar disk population, may show the opposite picture, with 
progressively lower eccentricities at larger distances from 
the MBH. For this reason the distribution of eccentric- 
ities vs. distance could be a key observational property 
for discriminating models for the origin of the young B- 
stars. We do note, however, that very close to the MBH 
(the innere regions of the S-cluster art a few 0.1 pes from 
the MBH), where general relativistic precession due to 
the MBH becomes important. At these regions resonant 
relaxation might be quenched, and capture stars could 
again show higher eccentricities. 

3.4. Mass function of the captured stars 

The mass function (MF) of captured stars is likely to be 
regular at the high mass regime (> 3Mq), i.e. reflecting 
the mass function of stars far from the MBH, where regu- 
lar star-formation can take place, and where the progen- 
itor binaries form. The currently observed MF, however, 
may slightly differ from the regular MF due to several 
reasons. 

(a) Binary fractions and the distribution of binary or- 
bital properties are mass dependent. Massive (OB) stars 
have higher binary fractions than lower mass stars. Mas- 
sive binaries are also typically much more compact (a 



log constant distribution of period s vs. a log-normal dis- 
tribu t ion for low mass stars; e.g . iDuquennov fc Mavod 
I199H iKobulnickv fc Frved l2007f i. The binary disrup- 
tion process requires binary progenitors for the captured 
stars. It also maps the binary period into the capture 
separation from the MBH as discussed above. For both 
these reasons, a larger fraction of massive stars is likely 
to be captured in this scenario, compared with the frac- 
tion of low mass stars binary disruption (relative frac- 
tions; in terms of absolute frequency of low mass binary 
stars, they are still more frequent than high mass binary 
stars). In addition, massive stars are typically captured 
much closer to the MBH than low mass stars. Therefore, 
the population of captured stars in the GC, especially in 
the central pc, tends to be biased towards massive OB 
stars. 

(b) Longer living stars (hence lower mass stars) cap- 
tured at earlier times may have a higher probability of be- 
ing dis rupted by t he MBH during their dynamical evolu- 
tion (Peret s et aLl[2009( l. On even longer timescales (two- 
body relaxation timescales), they can migrate from their 
original capture orbit. For this reason, the observed pop- 
ulation of captured stars is biased towards stars which 
were recently captured, i.e. selected against long living 
stars captured at earlier times. 

(c) A small contribution from Kozai induced merger 
of captured binary stars, induced by Kozai resonances 
near the MBH, m ay contribute a sm a ll fraction of 
more massive stars (jAntonini et aTl 120091 ; iPeretd l2009bl ; 
iPerets fc Fabrvckvll2009[ ). 

Currently, the MF of the probed B-stars population is 
consistent with a continuous star fo rmation (with regu - 
lar IMF) over the last ~ 50 Myrs (jBartko et al.ll2010f ). 
The range of probed masses is currently quite limited, 
and low mass stars are not included in the analysis (and 
therefore point (a) is not directly verified, given our lim- 
ited knowledge on the binary distribution of OB stars of 
different masses. For this limited range, the observed MF 
is qualitatively consistent with the expected one. Future 
observations probing fainter stars will provide much bet- 
ter data and therefore stronger constraints on the mass 
function. 

4. DISCUSSION 

Any model suggesting a common origin for all/most of 
the B-stars in the GC is required to produce both the 
stars in the inner region (< 0.05 pc) and those in the 
outer region (< 0.5 pc or more). Most of our discussion 
on the constraints on different models for the origin of 
these stars can therefore be summarized by two simple 
qualitative properties: (i) the ratio between the number 
of stars in the inner region and the outer region of the 
GC, (ii) the radial eccentricity distribution trend in these 
regions. 

Table [T] summarizes the predictions of the various the- 
oretical models for the number of B-stars in the in- 
ner/outer regions, compared to current observations. 
As can be seen, most models grossly over-predict the 
number of B-stars in the outer region, up to 0.5 pc, 
given the observed number of B-stars in the central 0.05 
pc. Such constraints could easily be studied and ex- 
tended/generalized for any new model suggested for the 
B-stars origin. 

Currently, we find that the most consistent explanation 
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Fig. 3. — The cumulative eccentricity distribution in three different distance regimes, after 20, 40 and 60 Myrs of evolution. Also showrn 
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population is ~ 0.7 (^/(l/2). Captured stars have initial median 
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for the complete population of B-stars in the GC is an 
external formation followed by binary disruption (likely 
driven by massive perturbers). The eccentricity distri- 
bution of the captured stars in this scenario should be 
dependent on the distance from the MBH. Higher eccen- 
tricities should be observed for stars captured at larger 
distances. 

The formation of binary B-stars in a very eccentric 
stellar disk with einit > 0.7 could also be consistent with 
the ratio of inner to outer B-stars, for a high enough bi- 
nary fraction in the stellar disk. The radial distribution 
of the stars in this case is likely to differ from that of the 
massive perturbers scenario, as the population of B-stars 
would include both the progenitor population formed in 
the stellar disk and the captured stars. These two pop- 
ulations are likely to have different radial distributions, 
and their combined distribution is likely to show a bi- 
modality, with one peak at the inner region close to the 
MBH, dominated by the captured B-stars, and an ad- 
ditional peak corresponding t o the peak in the ste llar 
disk. From the simulations bv iMadigan et ahl (|2009[ ). it 



seems that the binary disruption is only weakly depen- 
dent on the binary separation. This is reminiscent of 
the empty loss cone regime discussed above. The radial 
distribution of the captured stars should therefore cor- 
respond to a direct mapping of the binaries separations 
distribution in this case. The eccentricity distribution of 
the B-stars should also be a combination of two popula- 
tions. A more detailed study of this scenario, however, is 
beyond the scope of th is paper and will be studied else- 
where (see Mad igan et al.i i200S for a discussion of the 
basic scenario). 

5. SUMMARY 

In this study we explored the dynamical constraints 
on the origin of the young B-stars in the Galactic center. 
We showed that most of the currently studied models 
could not consistently explain both the population of the 
S-stars (the B-stars at < 0.05 pc from the MBH), and 
the extended population of young B-stars up to 0.5 pc. 
In these models the number of B-stars in the inner and 
outer regions are incompatible. 

We used analytical arguments and N-body simulations 
to study the capture and the dynamics of the B-stars 
in the external formation+binary disruption scenario, 
which is one of the models consistent with the observed 
B-stars numbers at different regions of the GC. Using our 
analysis we provide detailed predictions for the distri- 
bution of the B-stars dynamical properties, which could 
serve as observational signatures to further support, con- 
strain or refute this model. 
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